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We investigate the critical properties of two-dimensional Z(N) vector models for N larger than 4. 
In particular, critical points of the two phase transitions are located and some critical indices are 
determined. We study also the behavior of the helicity modulus and the dependence of the critical 
points on N. 
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1. Introduction 

The Berezinskii-Kosterlitz-Thouless (BKT) phase transition was originally discovered in the 
two-dimensional (2D) XY model in the first half of the seventies [1]. Since then it was realized that 
this type of phase transition takes place in a number of other models including discrete 2D Z(N) 
models for large enough N and even 3D gauge models at finite temperature (see [2, 3] for a recent 
study of the deconfinement transition in 3D t/(l) lattice gauge theory). Here we are interested in 
the phase structure of 2D Z(N) vector models. On a 2D lattice A = L 2 with linear extension L and 
periodic boundary conditions, the partition function of the model can be written as 



Z(A,j3) 



xeA ly s(x)=0_ 



exp 



2k 

L L P cos TF (s(x)-s(x+e n )) 



N 

xeA n=l,2 ly 



(1.1) 



The BKT transition is of infinite order and is characterized by the essential singularity, i.e. the 
exponential divergence of the correlation length. The low-temperature or BKT phase is a massless 
phase with a power-law decay of the two-point correlation function governed by a critical index rj . 
The Z(N) spin model in the Villain formulation has been studied analytically in Refs. [4]. It was 
shown that the model has at least two BKT-like phase transitions when N > 5. The critical index 77 
has been estimated both from the renormalization group (RG) approach of the Kosterlitz-Thouless 
type and from the weak-coupling series for the susceptibility. It turns out that T7 (j8i 1 ) = 1/4 at 

the transition point from the strong coupling (high-temperature) phase to the massless phase, i.e. 

(2) 

the behavior is similar to that of the XY model. At the transition point /3 C from the massless 
phase to the ordered low-temperature phase, one has 7](j3 c ) = 4/N 2 . A rigorous proof that the 
BKT phase transition does take place, and so that the massless phase exists, has been constructed 
in Ref. [5] for both Villain and standard formulations. Monte Carlo simulations of the standard 
version with N = 6,8, 12 were performed in Ref. [6]. Results for the critical index r\ agree well 
with the analytical predictions obtained from the Villain formulation of the model. 

In Refs. [7, 8] we have started a detailed numerical investigation of the BKT transition in 2D 
Z(N) models for N = 5 which is the lowest number where this transition can occur. Our findings 
support the scenario of two BKT transitions with conventional critical indices. Here we continue 
our study with investigation of models for N=7 and 17. We want to locate the transition points and 
to compute some critical indices. Such results could serve as checking point of universality in our 
studies of the BKT transitions in 3D gauge models at finite temperature. Our second goal is to use 
the available data on the position of critical points to deduce phenomenological scaling of these 
points with N. 



2. Numerical results 

We simulated the model defined by Eq. (1.1) using the same cluster Monte Carlo algorithm 
adopted in the Refs. [7, 8] for the case N = 5. We used several different observables to probe 
the two expected phase transitions. In order to detect the first transition (i.e. the one from the 
disordered to the massless phase) we used the absolute value \Mi | of the complex magnetization, 
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Figure 1: Susceptibility %i versus j3 in Z(7) (left) and Z(17) (right) on lattices with several values of L. 



and the helicity modulus [9, 10] 

Y = (e)-L 2 p(s 2 ), (2.2) 

where e = X H<ij>, cos (6/ — fy) an d s = ji H<ij> x sin(0,- — For the second transition (i.e. the 
one from the massless to the ordered phase) we adopted the real part of the "rotated" magnetization, 

M R = \M L \ca&(Ny) , 

and the order parameter 

niy = cos(N\j/) 

introduced in Ref. [11], where y is the phase of the complex magnetization defined in Eq. (2.1). 
In this work, both for N = 1 and N = 17, we collected typically 100k measurements for each value 
of the coupling /3, with 10 updating sweeps between each configuration. To ensure thermalization 
we discarded for each run the first 10k configurations. The jackknife method over bins at different 
blocking levels was used for the data analysis. 

In Fig. 1 we show the behavior of the susceptibility = L 2 {(\Ml\ 2 ) — (|M/J) 2 ) of the ab- 
solute value of the complex magnetization, which exhibits, for each volume considered, a clear 
peak signalling the first phase transition. The position of the peak in the thermodynamic limit 
defines the first critical coupling, J3 C . Fig. 2 shows instead the behavior of versus j8 on var- 

(2) 

ious lattice sizes; here the second critical coupling j8 c is identified by the crossing point (in the 
thermodynamic limit) of the curves formed by the data on different lattice sizes. 

To determine the first critical coupling J3 C , we could extrapolate to infinite volume the pseudo- 
critical couplings given by the position of the peaks of xf 1 " 1 . However, since the approach to the 
thermodynamic limit is rather slow (powers of logL), we adopted a different method, based on the 
use of the "reduced fourth-order" Binder cumulant 



3(|M L p 



U L - 1 0/1^12X2' 



the cumulant B^' defined as 



R (M R )_ (\Mr-(M r )\\ 
* 4 - (\M R -(M R )\^> 
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P P 
Figure 2: Behavior of m v with /3 in Z(7) (left) and Z(17) (right) on lattices with several values of L. 

and the helicity modulus T. We estimated j3c by looking for (i) the crossing point of the curves, 
obtained on different volumes, giving a Binder cumulant versus j8 and (ii) the optimal overlap of 
the same curves after plotting them versus (j8 — /3 f )(logL) 1//v , with v fixed at 1/2. The method (ii) 
has been applied also to the helicity modulus T. Our best values for jSj 1 ^ are 

N = 7 : j3 c (1) = 1.1113(13) , N = 11 : /3 f (1) = 1.11375(250) . 

Then, we performed the finite size scaling (FSS) analysis of the magnetization \Mi\ and the sus- 
ceptibility at Pc using the following laws: 

|M L |(/3 C (1) ) =AL-V V , *f >(/3 c {1) ) =BZ//v , (2.5) 

where y/v = 2 — 7] and tj is the magnetic critical index. Results are summarized in Tables 1 and 2. 
We observe that the hyperscaling relation y/v + 2/3/V = d = 2 is nicely satisfied within statistical 
eiTors in both models. 



Table 1: Results of the fit to the data of \M L \{& 1) ) with the scaling law (2.5-left) on L 2 lattices with L^L, 
for N = 7 andJV= 17. 



N = l 






A 


J3/V 


£ 2 /d.o.f. 


N=ll 


32 


1.00653(48) 


0.12210(08) 


5.5 


^min 


A 


J3/V 


X 2 /d.o.f. 


64 


1.00858(70) 


0.12243(12) 


3.7 


32 


1.00388(51) 


0.12111(09) 


7.98 


128 


1.01074(94) 


0.12277(15) 


2.0 
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1.00620(69) 


0.12149(12) 
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256 


1.0146(16) 


0.12336(26) 


0.40 
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1.0089(11) 


0.12191(18) 
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1.0162(22) 
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256 


1.0107(15) 


0.12219(24) 


0.74 


512 


1.0177(38) 


0.12381(56) 


0.16 


384 


1.0113(24) 


0.12228(36) 


1.36 
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1.0185(57) 


0.12393(84) 


0.28 





We can cross-check our determination of the critical exponent tj by an independent method, 
which does not rely on the prior knowledge of the critical coupling, but is based on the construction 
of a suitable universal quantity [12, 8]. The idea is to plot y^^L^ 2 versus B { 4 Mr) and to look for 
the value of tj which optimizes the overlap of curves from different volumes. We found that, both 
in Z(7) and Z(17), tj = 1/4 is this optimal value, since it gives the best overlap of these curves in 
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Table 2: Results of the fit to the data of Xl(P^) w ^ me scaling law (2.5-right) on L 2 lattices with 
L ^ L min , for N = 1 and N = 17. 



N = l 




^min 


B 


y/v 


% 2 /d.o.f. 




N 


= 17 




32 


0.00558(07) 


1.7402(23) 


2.38 


^min 


B 


y/v 


X 2 ld.o.f. 
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0.00540(09) 


1.7457(29) 


1.32 
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0.00559(08) 


1.7372(26) 


2.7 
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Figure 3: Correlation between %\ Mr) L^ 2 and the Binder cumulant B { ™ r) for 77 = 0.25 in Z(7) (left) and 
Z(17) (right) on lattices with L ranging from 128 to 1024. 

the region of values corresponding to the first phase transition, i.e. the lower branch of the curves 
of Fig. 3. This result for v\ agrees with the determinations v\ = 2 — y/v from the FSS analysis. 

As for the second critical coupling /3 C , we used the same method adopted for jSj 1 ', but applied 
now to B4 and m^. Our best estimates are 

Ar = 7: y3 t (2) = 1.8775(75) , N=\l: j3j 2) = 10.13(12) . 

The standard FSS analysis applied to the susceptibility of the rotated magnetization Mr at 

/3j 2 ^ leads to the result for the critical indices y/v given in Table 3 l . 

Also in this case the critical index 77 can be determined by an independent method, irrespec- 
tively of the knowledge of (5^: MrL^I 2 is plotted versus my and the value of 77 is searched for, 
which optimizes the overlap of data points coming from different volumes. The results we found 
for 77 in Z(7) and Z(17) are in perfect agreement with the theoretical prediction r\^> = 4/N 2 (see 
Fig. 4). 

Finally, in Fig. 5 we present the behavior with /3 of the helicity modulus (2.2). This quantity is 
constructed in such a way that it should exhibit a discontinuous jump (in the thermodynamic limit) 

1 We do not report in this work the determinations of fi/v by the FSS analysis of the rotated magnetization Mr, 
since they are affected by large statistical and systematic uncertainties. 
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Table 3: Results of the fit to the data of X^iPP) 
L ^ L mtn , for N = 7 and N = 17. 



with the scaling law (2. 5 -right) on L lattices with 
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Figure 4: Correlation between MrL^I 1 and m v for r\ = A/1 2 in Z(7) (left) and 77 = 4/17 2 in Z(17) (right), 
on lattices with L ranging from 128 to 1024. 
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Figure 5: Helicity modulus versus j3 in Z(5), Z(7) and Z(17) on lattices with various sizes. 



at the critical temperature separating the disordered phase from the massless one, if the transition is 
of infinite order (BKT). Since the Kosterlitz-Thouless RG equations for the XY model [1, 13, 14] 
lead to the prediction that the helicity modulus T jumps from the value 2/(71 j3 ) to zero at the critical 
temperature, one can check if the same occurs for vector Potts models. In Fig. 5 we plot a red line, 
representing the function 2/(71/3); the crossing between this line and the curves formed by data 
points of T approaches indeed jSj 1 ' when the lattice size increases. 

We studied also the specific heat at the two transitions, finding that, in contrast to the case of 
first- and second-order phase transitions, it does not reflect any nonanalytical critical properties at 
the critical temperatures, thus confirming that only BKT transitions are at work here. 
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3. Summary 

We have determined the two critical couplings of the 2D Z(N = 7,17) vector models and 
given estimates of the critical indices v\ at both transitions. Our findings support for all N > 5 the 
standard scenario of three phases: a disordered phase at high temperatures, a massless or BKT one 
at intermediate temperatures and an ordered phase, occurring at lower and lower temperatures as N 
increases. This matches perfectly with the N — > °° limit, i.e. the 2D XY model, where the ordered 
phase is absent or, equivalently, appears at j8 — > °°. Considering the determinations of the critical 
coupling j8 c ' (N) obtained in this work for N = 7 and 17, together with those obtained for N = 5 
in Refs. [7, 8] and those for ,/V = 6, 8 and 12 of Ref. [6], one can verify that ySj (N) approaches the 
2D XY value, j8 c (1) = 1 . 1 199, exponentially in N or even faster, while /3j 2 ^ (N) grows to infinity with 
N 2 . A more detailed analysis of the A^-behavior of critical couplings will be given elsewhere [15]. 

We have also found that the values of the critical index tj at the two transitions are compatible 
with the theoretical expectations. 
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